N91-14980 


STUDYING THE SPATIAL DISTRIBUTION OF INTERSTELLAR DUST 

H. Walker,*’** M.Werner,* C. Allen,* R. Henry,*** R.Kimble,*** J. Wofford,*** 
J.Murthy**** 


* NASA Ames Research Center, Moffett Field, CA 94035 
** SETI Institute, Moffett Field, CA 94035 

*** Center for Astrophysical Sciences, The Johns Hopkins University, Baltimore, 
MD 21218 

**** NASA Goddard Space Flight Center, Greenbelt, MD 20771 


I. Introduction 

The spatial distribution of interstellar dust reflects both interstellar dynamics and 
the processes which form and destroy dust in the interstellar medium. The IRAS 
survey, because of its high sensitivity to thermal emission from dust in the infrared, 
provides new approaches to determining the spatial distribution of the dust. We 
report here the initial results of an attempt to use the IRAS data to probe the 
spatial distribution of dust - by searching for thermal emission from dust in the 
vicinity of bright stars. 

It can easily be shown that a bright, luminous star embedded in a region of “average” 
interstellar density, i.e. n(H) = 1, and “typical” dust-to-gas ratio, will heat the 
nearby dust and produce infrared emission readily detectable from IRAS. Consider 
for example a B3 star at a distance of 50pc from the earth. At a radial distance of 
0.25 pc from the star (angular distance 0.25° as seen from earth), the heat input 
into the interstellar dust from this star is 1.1 x 10 -21 erg cm 3 s 1 H-atom 1 , more 
than 200 times that due to the “mean” interstellar radiation field. Assuming n(H) 
= 1 and a path length of 0.25 pc through the material within 0.25 pc of the star, the 
reradiated infrared luminosity due to this heat input will be 6.6 x 10 5 erg cm 2 sr 1 
s -1 . If 50% of this reradiated power lies in either of the two long wavelength IRAS 
bands, the resulting surface brightness will be l-3MJy/sr, detectable by IRAS in 
regions where the background is smooth, due to the large number of pixels involved. 
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II. Results to Date 

The above analysis has motivated us to search for diffuse infrared emission in the 
vicinity of bright, luminous stars, in order to assess the probability with which 
they are associated with interstellar dust clouds in certain density ranges. We 
began with a list of 25 stars selected on the basis of high visible brightness or high 
ultraviolet brightness, with |b| >15°. The 60^m and lOO/xm IRAS Sky Flux plates 
from HCONS 1 & 2 (resolution 6 arcmin) were examined in a search for infrared 
emission associated with these stars. In five of the twenty-five cases, extended 
infrared emission was found close to or at the position of the star. In each case, 
this emission shows 60pm/100/xm flux ratio appreciably higher than that which 
characterizes the diffuse emission elsewhere on the same plate. The higher color 
temperature indicates that the infrared emission is due to dust heated by the star 
rather than background or foreground “cirrus” clouds. 

III. Analysis 

We list below the properties of the diffuse infrared emission associated with these 
five stars, a Lyr (Vega) is included in the table as a comparison object, known 
to have a very cool circumstellar dust shell, and one of the IRAS standard point 
sources, a Boo (Arcturus), is also listed. The table gives the peak surface brightness 
(Ipea*) of the emission at the location of the star, with the background subtracted 
(in MJy/sr). The background is found by taking a ring around the area used for the 
summation and finding the mean of the lowest 25% of the values. This avoids the 
problems of nearby sources of emission and of stripes, and since the areas involved 
are small, the assumption of a single value for the background is acceptable. The flux 
density, with the background subtracted, is summed over a region usually 22 arcmin 
by 22 arcmin (X) I net in Jy). The values are shown in brackets when the emission 
peak is not precisely centred on the source position. For a Vir two boxes were used 
to sum the emission, one centred on the 12/zm source and one on the lOO/im source. 
One box was used for 6 Sco, but its size was increased to 48 arcmin by 38 arcmin to 
accomodate the extended emission at lOO^m. The 60^m/l00/zm flux density ratio 
is shown in two forms. The subscript “net” means that the preceeding two rows of 
summed background subtracted flux are divided to give the ratio. The subscript 
“sub” means that the two maps with the backgound subtracted from each were 
divided and the value of an appropriate pixel given. The same ratio for cirrus is 0.2 
(Weiland et al., 1986). The numbers in brackets are the temperatures derived from 
the ratio, assuming a blackbody. The ratio of the infrared luminosity to the stellar 
luminosity at earth is given, showing how small an optical depth is detected. 

Also given is the heating rate 1 parsec from the star (r). The stellar input spectra 
are derived from Kurucz (1979) models (a blackbody is used for a Boo). These 
spectra are then folded through the absorption properties of dust in the model 
by Mathis, Mezger and Panagia (1983) to determine the heat input, assuming the 
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dust is optically thin. The canonical relationship between extinction and hydrogen 
column density is assumed, E b-v = Njy/6xl0 . 

IY. Conclusions 

These preliminary results show that this technique - which relies on finding infrared 
emission associated with randomly selected stars - can ultimately be used to study 
the distribution of dust in the ISM. 

The density of the cloud producing the infrared emission may be derived by assum- 
ing that the dust is at its projected distance from the star and that the heating 
is due to the star’s (known) radiation field. The heating radiation will be folded 
into a Draine (1985) grain model, and the number of emitting grains adjusted to 
reproduce the observed energy distribution. We note that this technique is capable 
in principle of detecting dust densities much lower than those typical of the cirrus 
clouds, because we are looking preferentially at regions near stars where the heating 
flux is far higher than the diffuse radiation field which produces the cirrus emission. 
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Notes for the Table 

a Vir is also classified as B1III-IV. The fluxes for a Vir in brackets show that the 
flux given is not centred on the star, but about 14 arcmin away. The 12^/.m/25/^m 
map shows a single peak at the stellar position, showing the star dominates the 
emission here. The 60//m/l00^m shows two features, one at the stellar position 
and one at the peak of the 100/xm flux. The flux ratios include both the star and 
the dust. 

For 6 Sco the peak 12fim flux from the star is not coincident with the peak of the 
100/im flux map. Van Buren and McCray (1988) suggest the longer wavelength 
emission is a bow shock associated with the star. In the 12^xm/25p.m map the star 
is at a peak in the ratio map, and the cool dust is near a minimum. The structure is 
poorly defined in the contour plot for the 60/im/l00/zm map, but the bow shock 
structure is present. 

v Sco has a companion about 40 arcsec away, which is itself double. The 12pm/25/mi 
map shows a minimum in the ratio at the star’s position, and there is a maximum 
in the 60^m/100//m map, as the cool dust dominates the surroundings. 

7 r Sco is visible at all 4 wavelengths, and appears much larger than a point source. 
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The 12/im/25//m shows a minimum in the ratio at the source position, and the 
60^m/l00^m map shows a peak. 

The luminosities of the stars, shown in the table, are from Allen’s Astrophysical 
Quantities. 

r is the heating rate 1 parsec from the star. 
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